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The Jarzynski equality (JE) is known as an exact identity for nonequillib- 
rium systems pp. The JE was originally formulated for isolated and isothermal 
systems, while Adib reported an JE extended to an isoenergetic process[2]. 
On an isoenergetic process, the functional form of the Hamiltonian is altered 
with the energy fixed at a constant value E under the effect of an aritificial 
field term which is added to the Hamilton dynamics. This kind of JE is useful 
to calculate the entropy which replaces the Helmholtz free energy appearing 
in the original JE. 

On the other hand, the limitation of the constant energy seems to be an 
obstacle to widely apply the JE to more interdisciplinary problems such as 
optimization problems. Another problem is that it takes too long time to 
estimate some physical quantities since the entropy at the different energy 
values cannot be calculated at the same time with the JE. In the present 
paper, we extend the JE to an energy-controlled system in order to overcome 
these difficulties. In our study, we make it possible to control the instanta- 
neous value of the energy arbitrarily in a nonequilibrium process. Under our 
extension, the new JE is more practical and useful to calculate the number 
of states and the entropy than the isoenergetic one. Furthermore, we expect 
that our JE can be used also for applications to computation and information 
science beyond physics. We will show an instance of possible application of 
our JE to a kind of optimization problems which is equivalent to physical 
problems to find the ground state [3]. 

Following the derivation of Adib, let us consider a classical system which 
evolves under the following dynamics, where an artificial field term F = 
{Fx, Fp) is added to the Hamilton dynamics. 
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The above T = {x,p) describes a point on the phase space. In order to 
control the value of the energy as a function of time E{t) from t = to r, 
we choose the functional form of -F = {Fx, Fp) as 

where X is a vector on the phase space satisfying X-dvH ^ 0. From eqs. ([I]) 
and (|2]), we immediately recognize that the energy is accurately controlled 
as E{t) since we have 

dH _dH ■ OH _ dE 

~dt~w^~dr~^' ^ ^ 

Under the dynamics ([1]), the ensemble density ptiX) evolves following the 
Liouville equation 

p,(ro = po(ro)e-*^(^'), (4) 



where 
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AtiTt) = - I dt'A{Tt,). (5) 







Equation (|5]) is the time average of the "phase space compression factor" 
A(rt) = drt ■ Tt along the trajectory that connects Fq to rt.[4j This factor 
determines the time evolution of the ensemble density pt(rj) from t = to 
r, and plays the central role to estimate the entropy efficiently. The system 
is set to be in an equilibrium state at the initial time t = similarly to the 
original JE. Therefore po(ro) is equal to the microcanonical distribution at 
E = E{0) 

"0 

where Qq is the number of states at t = 0. 

Let us consider the average of e"^^^ over all possible realizations from t = 
to t: 

(e'^) = j dr,p,(r,)e'^ (7) 

dr/'^'^°>-^'°)>. (8) 
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In the second line, eqs. (jl]) and (jS]) have been used. Since the value of the 
Hamiltonian is controlled as H{Ft) = E{t), we have 

H{rr)-E{r) = H{ro)-E{0). (9) 

Therefore -f^(ro) —E{0) in the Dirac delta function of eq. dH]) can be replaced 
by H{Tr) — E{t). Then, eq. ([S]) is reduced to 

(e-^) = (10) 

where AS" = ln(fi^/f2o) is the entropy difference between the equilibrium 
states at the different energy values E{0) and E{t). If we estimate the en- 
tropy difference with such a naive method as directly calculating the entropy 
at the initial and final times, we have to make the estimation of entropy 
twice. On the other hand, our formula ( ITOl) enables us to obtain the entropy 
difference by taking just a single average of e"^^^ . The above equality is the 
our central result, and we will show a practical aspect of our formulation 
next. 

In the application to an optimization problem, let us consider a simple 
optimization problem to find a minimum of the potential energy. If the 
potential energy has a complex structure, simple searching techniques such 
as the steepest decent method[5] does not work well to reach the minimum[6]. 
Therefore we need to employ an approximate algorithm in practice. We show 
that our JE can be used for a quantitative estimation, which indicates how 
much a tentative solution given by the approximate approach differs from 
the actual minimum solution. 

Let us consider an arbitrary form of the potential energy V{x) of contin- 
uous variables x which has no degeneracy at the ground state, and assume 
that an approximate solution with energy £ has been obtained through an 
approximatiove algorithm. With the theory to control the energy, the energy 
of the system is changed from -E'(O) = S' to E{t) = £, where £' is larger than 
S. As the Hamiltonian H is set as 

H = p' + -V{x), (11) 

r 

we can easily calculate the entropy at the initial time ^(O) since the system 
equals to non-interacting particles with the mass of each particles m = 1/2 
at t = 0. Then, the entropy at the final time S{t) = S{0) + AS is obtained 
with AS estimated in the following way, and the results are independent of 
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Figure 1: The energy is controlled from £' to £ for the potential energy 
V{x) = — x/2 — 2sin(37r(x + 2)) — 1/2. Length of the horizontal lines 
indicating £ and £' corresponds to the entropy at each energy. 

the value of r. First, we randomly choose an initial condition {xq,Pq) from 
the set {(a;, p)\H{x,p;t = 0) = £'} and note that the initial Hamiltonian 
depends only on p. We obtain the path toward a phase point with the 
lower energy following the energy-controlled dynamics ([1]) under the given 
initial condition. Some initial conditions result in the divergence of the factor 
appearing on the right-hand side of eq. (|2]) when the system is trapped in a 
local minimum with drtH = 0. Such a divergence means that the energy is 
not to be decreased to £ and the system evolving from such initial conditions 
are unable to reach any points on the phase space at t = t. Therefore such 
samples should be excluded in taking the average of e'^^^. Finally, we take 
the average of e"^^^ only in the case of the absence of such divergences, and 
we obtain AS appearing on the right-hand side of our JE ffTOj) . 

The minimum value of entropy is —oo in the case of classical systems. 
Therefore S{t) — —oo implies that the solutions for the given potential 
energy are close to the minimum point. In other words, estimating the dif- 
ference of the entropy can be an indicator how much the solution differs from 
the actual minimum. 

For instance, we estimate the entropy gap AS* for the potential 

V{x) = - x/2 - 2 sin(37r(x + 2)) - 1/2, (12) 

where r = 10, £' = 3, and £ = 2.9, 2.99 as Fig . 1. Results are AS = 0.561 
{£ = 2.9), 0.464 = 2.99). A lower £ = 2.9 gave a larger entropy difference 
0.561. The difference of AS between 0.561 and 0.464 is equal to 0.097, which 
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indicates that the case oi S = 2.9 is closer to the actual minimum solution 
than the case oi £ = 2.99 by 0.097. 

JE f lTOj) can be used also for the calculation of the canonical average of 
an observable 0{E), which can be expressed as follows, 

(^(^)) - JdEe3(^)~P^ ' ^^^^ 

where e'^*^^-' is equal to the number of states at the energy E and can be 
obtained from the JE (fTOj) . While Adib has introduced a similar method 
with JE for an isoenergetic process[2j, such a method takes a long time in 
practice to estimate the entropy as a function of E. In our method we can 
control the energy dynamically and therefore obtain the values of e'^'^^ at 
various values of the energy from E{0) to E{t) in a single realization. Then, 
the necessary time to take the above average can be much shorter than the 
method of Adib. 

In future, we would modify the JE on an energy-controlled process for a 
quantum system. A quantum system is more useful to estimate the efficiency 
of an approximate solution of an optimization problem since the entropy at 
the ground state is equal to which is clearly more suitable for the quanti- 
tative estimation than — oo of classical systems. Furthermore, the tunneling 
effect and discrete energy levels seem useful also for estimating the entropy 
of systems with discrete variables or complex landscape of potential energy 
like spin systems. 
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